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ABSTRACT 

We present faint galaxy counts from deep VRI images obtained with the Keck Telescope. 
These images reach R ~ 27 in median seeing FWHM ~ 0.5-0.6 arcsec and we detect a 
integrated galaxy number density of 7 x 10 5 degree -2 , equivalent to 3 x 10 10 galaxies in the 
observable Universe. In addition we present median galaxy colors as a function of magnitude; 
bluing trends are visible in all colors to R ~ 24.5. Fainter than R ~ 24.5, however, the typical 
(V—R) color becomes redder again, (V — I) remains constant, and (R—I) becomes yet bluer. 
These trends are consistent with the VRI count slopes, implying a decrease in the V slope 
at the faintest levels, which our data supports. Taking advantage of our good seeing we 
also present median half-light radii for faint galaxies, these show a steady decline at fainter 
magnitudes, leading to an intrinsic half-light radius of ~ 0.2 arcsec for a typical R ~ 25.5-26 
galaxy. Irrespective of the redshift distribution, the extremely high galaxy surface densities 
and their small intrinsic sizes are consistent with a scenario in which the majority of the very 
faint field population are dwarf galaxies or sub-galactic units. 

Subject headings: cosmology: observations - galaxies: evolution - galaxies: photometry 
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Introduction 



The issue of the total number of galaxies in the visible Universe is more than an academic curiosity; 
number counts as a function of apparent magnitude hold important clues to the formation and evolution 



of galaxies (see Koo & Kron 1992 for a review). The study of very faint field galaxies using CCDs was 
pioneered by Tyson (1988) who found an integrated number density of ~ 2 x 10 5 degree -2 , or ~ 10 10 
galaxies over the whole sky at B ~ 26. Two more recent studies which are notable for their depth and 



resolution are those of the Hawaii group (Lilly et al. 1991) and the Durham group (Metcalfe et al. 1991 



Metcalfe ct al. 1995|) . The latter now reach galaxy surface densities in excess of 4 x 10 5 degree 2 for a 



small field at a limiting magnitude of B — 27.5. Equivalent surface densities are also being achieved in the 
near-infrared, K ( pardner et al. 1994 , Cowie et al. 1995, Djorgovski et al. 1995). 



Potentially the most interesting result from these studies, first observed by Lilly et al. (1991) and 
confirmed by Metcalfe et al. (1995) with deeper data, is the change in the slope, a (where log dN/dm ~ a m), 



1 Based on observations obtained at the W.M. Keck Observatory, which is operated jointly by the California 
Institute of Technology and the University of California. 



-2 - 



of the deep B counts. The differential galaxy number counts in the B-band show a decrease in slope from 
as ~ 0.5 to ~ 0.3 at B ~ 25 (equivalent to i? ~ 24-24.5). In addition, Roche et al. (1993) claim 
that the amplitude of the two-point angular correlation function reaches a minimum around B ~ 25. 
Taken together these two observations have been interpreted by Metcalfe et al. (1995) as a signature of an 
increasing dominance of low luminosity galaxies at z > 1. These features, if confirmed and observed in other 
passbands, will constrain galaxy evolution models, although redshift information will be required before 
any strong conclusions can be drawn. 

This letter presents the first deep optical counts with the 10-m Keck Telescope. These VRI 
observations are both deeper, at a given completeness limit, than previous published observations and more 
significantly were obtained in good seeing, FWHM ~ 0.5-0.6 arcsec, providing a wealth of information on 
the very faint field population. 



2. Observations and Analysis 



2.1. Observations and Reduction 



The observations presented here were obtained by S.R. Kulkarni and J.G. Cohen to study the optical 
emission from two high galactic latitude pulsars. Thus for our purposes the fields provide randomly selected 
samples of the faint field population. These data were obtained on the night of 1994 August 8 using the Low 
Resolution Imaging Spectrograph (LRIS, pke et al. 1995] ) on the 10-m Keck Telescope, Mauna Kea. The 
imaging mode of LRIS provides a 5.7 x 7.3 arcmin field onto a thinned Tek detector with 0.21 arcsec/pixel 
sampling and 8.0e~ readnoise. Conditions during this run were good and the night was photometric; a log 
of the observations is given in Table 1. The raw frames, with individual exposure times between 0.5-1.5ks, 
were dcbiascd and then flatfielded in a standard manner using dome flatficlds. For each passband the 
dithered frames, with pointing offsets of 10-20 arcsec, were registered and combined using standard IRAF 
routines and a clipped-average algorithm to produce the final images. We reproduce a randomly chosen 
lxl arcmin region from the R frame of 1640+22 in Figure 1. 

Photometric calibrations were derived from observations of Landolt (1992) standard stars at air masses 
very close to those of our science observations. We estimate the average errors in our absolute calibrations 
to be S < 0.02. Color differences between the standard stars and typical faint galaxies give systematic errors 
of a similar size. The estimated reddening is E(B — V) = 0.07 ( Stark et al. 1992] ) for both fields and so we 



have applied zero point corrections of SV = 0.16, SR = 0.12 and 51 = 0.07. We conservatively estimate our 
absolute magnitude calibrations to be good to <; 0.1 mag. 



2.2. Analysis and Photometry 



To analyse these images we elected to use the SExtractor image analysis package ( Bertin 1994 ). 
This package is a fast, image analysis program which can detect (using a standard isophotal limit and 
area cut algorithm), robustly deblend and analyse sparse and moderately crowded galaxy fields. The 
program provides positions, shapes, image profiles, half-light radii, isophotal and aperture magnitudes for 
all deblended objects detected on a frame. To detect objects we first smooth the frame with the point 
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spread function, then thrcsholded it at a level of 2.5a of the sky noise in the raw frame and identify objects. 
This threshold provides a typical point source detection limit of R ~ 27.2 and appeared optimal in our 
incompleteness and noise object analyses. 

To obtain total magnitudes for our galaxies we have followed the approach used by Djorgovski et 
al. (1995). For the fainter objects we use a fixed aperture (1.5 arcscc diameter or ~ 2-3xFWHM) to 
measure magnitudes, these are then aperture corrected to a 5.6 arcsec diameter under the assumption 
that the objects have roughly stellar profiles, a correction of S ~ —0.14 at R = 26. This assumption was 
tested by stacking many faint objects to make "average" faint galaxies with high enough signal-to-noise 
for large aperture photometry; we find that the average galaxies have aperture corrections close to those 
of the stars (as expected from their typical sizes, see below). For the brighter objects in the field, those 
with isophotal diameters larger than our adopted fixed aperture, we measure isophotal magnitudes above a 
surface brightness threshold of 0.9a of the sky noise (Table 1) and similarly correct to 5.6 arcsec diameter. 
The majority of the objects with magnitudes brighter than R ~ 20.5 are saturated in our frames; we do 
not present counts brighter than this limit. To measure colors for our objects we smooth all the frames to 
the effective seeing of the worst frame (the V data), before measuring photometry in 1.5 arcsec diameter 
apertures for all the objects selected from the R catalogs. 

Removal of stars from our object catalog is necessary at brighter magnitudes. With the good seeing 
experienced during these observations it is relatively straightforward to use image concentration to identify 
the stellar locus to a depth of R = 24.5. Comparison of the colors of these candidate stars to those in 
Landolt (1992) confirms their identification and supports our estimates of the errors on our absolute color 
calibrations of £ 0.1 mag. At R — 24.5 the stellar fraction is & 6% and falling, so we do not apply any 
correction for stellar contamination fainter than R — 25. 



2.3. Completeness Modelling 

Non-detection, failure to detect an object that is really there; false detection, detection of non-real, 
noise objects; and magnitude errors, incorrect measurement of an object's flux, all plague faint galaxy 
counts. These problems become progressively worse at fainter magnitudes, with their relative importance 
depending upon the slope of the counts (because magnitude errors "scatter" plentiful faint objects into 
brighter magnitude bins; Eddington bias) and the details of the data reduction procedures. 

We assess non-detection and magnitude errors simultaneously with a Monte Carlo simulation that 
involves adding artificial galaxies to the data and re-applying the detection algorithm. In order to simulate 
both the mean properties and morphological diversity of faint galaxies we create artificial objects by 
extracting galaxy images at magnitude levels for which errors are negligible and dimming them by 2.5 
magnitudes, at a constant angular size. It should be noted that if the scale sizes of galaxies decrease at 
fainter magnitudes (c.f. §3) then this procedure will tend to overestimate our incompleteness. Extensive 
simulations (> 10 objects per frame, dropped in one at a time) were used to create a large matrix Pjj 
each element of which is the probability that a galaxy with magnitude rrij is detected with a measured 
magnitude m^. This matrix contains all of the completeness and magnitude error information. On the 
assumption that the counts do not turn over quickly (i.e. a does not change rapidly with magnitude near 
the limit of the data), Py can be projected into a total detection rate rj, the ratio of number observed to 
true number at magnitude m, . The true number of galaxies in bin i can then be estimated by dividing the 
observed number Ni by the detection rate rj. 

To correct for false detections we have created noise frames identical to our observations, geometrically 
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remapped, combined and analysed them in the same manner as the observations. We subtract the number 
counts of noise objects found in these simulated frames from our observed counts. This correction is a few 
percent or less, even in our faintest magnitude bins. 

The faintest magnitude at which the counts can be reasonably completeness-corrected is subjective; 
we truncate our counts at or before the 50% completeness level, beyond which the Poisson errors, or any 
other error estimates, lose their meaning. 



3. Results and Discussion 



We plot in Figure 2 the raw and corrected differential galaxy number counts from our two fields. We 
measure slopes from the combined field counts of: <xy — 0.404 ± 0.015 for V = 22-24.25, flattening to 
a v = 0.28 ± 0.05 at fainter magnitudes, a R = 0.321 ± 0.001 for R = 21-27 and aj = 0.271 ± 0.009 for 
I = 19.5-25.5. We see no obvious variation in the R and I slopes. The evidence for the break in V count 
slope at V ~ 24.5 is certainly not pronounced, but it is stronger in the deeper, wider-area 1640+22 data, 
and as we show below it is also seen in the median galaxy colors. The integrated, corrected number counts 
are 2.9 x 10 5 per sq. degree to V = 26.5, 7.3 x 10 5 to R = 27, and 4.1 x 10 5 to I = 25.5. We also plot on 
Figure 2 counts by other workers. The discrepencies in the count slopes at bright magnitudes may arise 
from other groups choosing fields devoid of bright galaxies, this biases them towards underdense regions 
producing a systematic undercount of galaxies at bright and intermediate magnitudes and thus a steeper 
slope to the counts. At faint limits we find reasonable agreement between the slopes of the various datasets, 
with residual differences in normalisation consistent with the different photometric corrections applied by 
different workers and our quoted photometric errors. 

The color distributions for our faint galaxy sample are shown in Figure 3. All show a bluing trend with 
fainter magnitude until R ~ 24-24.5, where a typical galaxy has a flat spectral energy distribution in v f v 
(corresponding to (V — R) ~ 0.5 and (R — I) ~ 0.5). Fainter than R ~ 24-24.5 the typical (V — R) colors 
turn-around and become redder, while the (V — I) colors flatten out and the (R—I) continue to become 
bluer (c.f. |Steidel fc Hamilton 1993|) . If we are observing the same population in all three bands, then these 



trends in median color should be reflected in different count slopes in VRI. At R <, 24 the bluing trends 
indicate ay > O-r > oti as observed. After the turn-around, the colors imply a decrease in ay from steeper 
than / at the bright end, to matching the / slope at the faint end. This change in F-band slope is indeed 
seen in the counts, although its significance is not high. 

The good seeing of our observations allow us to study the angular sizes of the faint galaxy population. 
We plot in Figure 4 the observed half-light radii, rtu, measured from the i?-band images. A gradual, near 
linear, decline is seen in the median size of faint galaxies to R ~ 25.5-26, at which point the median size is 
only just distinguishable from the stellar locus. Representing the faint field population as exponential disks 
leads to an estimate of the intrinsic half-light radius of a typical R = 25.5-26 galaxy of rni ~ 0.2 arcsec. 
This angular size corresponds to less than 2.9 /i^q 1 kpc at any redshift in any world model, so it is clear 
that the typical objects we are detecting are either intrinsically smaller or produce a larger fraction of their 
emission in their nuclear regions than local bright galaxies. Splitting the sample in two in each magnitude 
bin on the basis of (V—I) color shows no statistical difference between the sizes of the blue and red galaxies. 

In conclusion, we have presented deep galaxy counts in VRI passbands. We find a decrease in slope in 
our bluest passband, V, at the faint end similar to that seen in the B-band by Metcalfe et al. (1995) at an 
equivalent apparent magnitude. Beyond R ~ 24 it appears that galaxy counts approach logdN/dm ~ 0.3m 
in all bands, shallow enough to remove the threat of a divergence in the integrated night sky brightness. If 
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we are not seeing to very high redshift, either because the Universe is Einstein-de Sitter and we are running 
out of volume, or because there are few luminous, high redshift field galaxies, then the common faint-end 
slope in the galaxy counts is a measure of the shape of the galaxy (or possibly pre-galactic fragment, in view 
of the small angular size of the faintest objects) luminosity function at the faint end. In these scenarios, the 
slope of 0.3 implies a faint end luminosity function Q(L) ~ L~ 1,75 at large look-back times, in contrast to 
L^ 1 observed locally (e.g. Loveday et al. 1992). The faint counts may thus confirm, at higher redshift, the 
steepening of the faint-end slope of the luminosity function suggested by analyses of spectroscopic samples 
at brighter magnitudes (Eales 1993, Ellis et al. 1995). Interestingly, the roll-over at V ~ 24.5 is within reach 
of the new generation of large telescopes and therefore we may hope for direct spectroscopic observations of 
galaxies in this intriguing magnitude range. 
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Tables 

Table 1 The table gives the log of the observations including the field identification and coordinates. For 
each filter we also quote the total exposure time, T exp ; the total number of sky photons detected per pixel, 
iV 7 ; the la surface brightness limit, /i(ler); the 50% completeness limit from our simulations, mu m ; the 
number of objects detected above this limit on the frame, Nu m , the FWHM of the seeing in arcsec and the 
total field area in sq. arcmin. 

Figures 

Figure 1 A randomly selected lxl arcmin field from the R image of 1640+22. Objects with magnitudes 
in the range R = 26.0-26.5 are marked. Notice the minimal crowding in the frame. 

Figure 2 Plots of differential galaxy counts as a function of magnitude in the VRI passbands in the two 
fields. The raw and corrected counts are represented with skeleton and solid symbols respectively, circles 
for field 1640+22, squares for 2229+26. The error bars are Poisson plus (on the solid points) an estimated 
uncertainty in the completeness correction. Least-square linear fits are shown with solid lines (there are two 
fits to the V counts split at V = 24.25). Where applicable, the work of other authors is shown with open 
points, these are truncated at the same completeness limit as our data. The following corrections have been 
applied to transform the different passbands onto our system: Lilly et al. (1991) I = Iab — 0.48 and Steidel 
& Hamilton (1993) R = K A b + 0.16. 

Figure 3 The typical colors of the faint galaxy population as a function of R magnitude in the two fields. 
The points arc medians of samples of 401 galaxies in each bin. The horizontal error-bars show the extent 
of the magnitude bin, while the vertical error-bars are 95% confidence limits calculated using boot-strap 
resampling of the data, non-detections are included, la limits within our photometry aperture lie outside 
the region of the color planes plotted. Note the turn-around in the (V — R) colors at R ~ 24.0-24.5 and the 
flattening of the (V — I) colors at the same magnitude. 

Figure 4 The half-light radii as a function of apparent R magnitude for all the objects in our two fields. 
The points are medians of samples of 401 galaxies with error-bars defined in a similar manner to Figure 3. 
The seeing difference between the fields has been removed by shifting the 2229+26 points. The dotted line 
marks the locus of stars, these are excluded from the medians. Each of the dashed lines is a linear fit to the 
variation of apparent half-light radius with source magnitude for an exponential disk with a fixed, intrinsic 
half-light radius, labelled in arcsec. The solid line shows the detection limit as a function of half-light radius 
for an exponential disk. For the most extended sources (rhi ^ 0.6 arcsec) a bias is visible, such that the 
source sizes are progressively underestimated at fainter magnitudes. This bias arises from a combination of 
incompleteness beyond R ~ 25 and the use of inappropriate aperture corrections, the corrections applied 
are calculated for "typical" , hence compact, faint galaxies. Nevertheless, it is apparent that the median 
source size in the field shrinks more rapidly than expected for a population with a fixed intrinsic angular 
size. By R ~ 25.5-26 the half-light radii is asymptotically approaching the stellar locus and indicates a very 
small intrinsic size for these faint galaxies, rhi ~ 0.2 arcsec. 
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